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The heterotrimeric enzyme AMP-activated protein kinase (AMPK) is a major metabolic factor that 
regulates the homeostasis of cellular energy. In particular, AMPK mediates the insulin resistance that is 
associated with type 2 diabetes. Generally, cellular processes require tight regulation of protein kinases, 
which is effected through their formation of complex with other proteins and substrates. Despite their 
critical function in regulation and pathogenesis, there are limited data on the interaction of protein kinases. 
To identify proteins that interact with AMPK, we performed large-scale affinity purification (AP)-mass 
spectrometry (MS) of the AMPK-al and -01 subunits. Through a comprehensive analysis, using a 
combination of immunoprecipitaion and ion trap mass spectrometry, we identified 381 unique proteins in 
the AMPKa/p interactomes: 325 partners of AMPK-al and 243 for AMPK-pl. Further, we identified 196 
novel protein-protein interactions with AMPK-al and AMPK-pi. Notably, in our bioinformatics analysis, 
the novel interaction partners mediated functions that are related to the regulation of actin organization. 
Specifically, several such proteins were linked to pancreatic beta-cell functions, including 
glucose-stimulated insulin secretion, beta-cell development, beta-cell differentiation, and cell-cell 
communication. 



Mammalian adenosine monophosphate-activated protein kinase (AMPK) is a serine-threonine kinase 
that regulates energy homeostasis at the cellular and whole-body levels 13 . AMPK is a heterotrimer that 
comprises a catalytic oc subunit, a scaffolding P subunit, and a regulatory y subunit 1,2 . Orthologs (otl, a2, 
pi, P2, yl> y2, y3) of AMPK subunits are found in all eukaryotic species and expressed in many tissues 2 , 
suggesting that the structure and regulation of AMPK are evolutionarily conserved 4 . 

AMPK is activated by nutrient deficiency and cellular stresses, such as glucose deprivation and hypoxia, which 
initiate metabolic adaptation programs to preserve cellular energy and maintain viability in various tissues, 
including skeletal muscle, cardiac muscle, adipose tissue, pancreas, and liver 2,3,5,6 . In general, AMPK activity is 
regulated through allosteric activation by adenine nucleotides, phosphorylation via upstream kinases 7,8 , and 
inhibition of dephosphorylation by protein phosphatases 9 . 

After activation, AMPK phosphorylates many substrates that regulate metabolism 10 , signaling 11 , gene express- 
ion 12 , and tissue-specific glucose and lipid homeostasis 13 . Based on these functions, AMPK is linked to various 
metabolic diseases, such as type 2 diabetes, obesity, hormonal disorders, cardiovascular disease, and cancer 2,7 . 
Drugs that activate AMPK, such as metformin and thiazolidinediones, are mainstays in the treatment of type 2 
diabetes, underscoring the importance of targeting the AMPK pathway to control metabolic disorders 2,3,5 " 8 . 

Recent evidence suggests that AMPK mediates many functions in pancreatic beta-cells, including glucose- 
stimulated insulin secretion (GSIS) 14 , proliferation 15 , and survival 16 . Because pancreatic beta-cells have distinct 
processes with regard to glucose metabolism and depend heavily on glucose and energy-sensing for their function 
and survival, the functions of AMPK must be examined in beta-cells, which have been described in liver and 
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skeletal muscle. Further, recent findings have demonstrated that 
AMPK is a positive and negative regulator of insulin secretion 14 . 
Based on this duality, considerable effort is under way to determine 
the functions of AMPK in pancreatic beta-cells. Yet, how AMPK 
regulates beta-cell function remains unknown 14 . 

Protein-protein interactions are crucial in all cellular processes 1718 . 
As expected, interactions between AMPK and proteins in its down- 
stream pathway are important for most of its functions 19 . Despite the 
significance of AMPK in beta-cell functions, the interactions between 
AMPK and other proteins in beta-cells have not been studied exten- 
sively. Thus, identifying target molecules that interact with AMPK is 
critical to understand the functions that are mediated by AMPK. 
Also, such information can increase our understanding of disease 
mechanisms and provide the basis for more specific therapies. 

Recent protein interaction studies that are based on mass spectro- 
metry-based proteomics have examined this issue in a systematic and 
relatively unbiased manner 17,18 , based on technological advances in 
mass spectrometry. Specifically, affinity purification, followed by 
mass spectrometry analysis (AP-MS), has become the preferred 
method of characterizing protein interactions 20,21 . 

We performed an extensive study of proteins that interact with 
AMPK-al and -(31, which are more highly expressed in pancreatic 
beta-cells than AMPK-a2 and |32 15,22 . To identify proteins that inter- 
act with these subunits, we performed an affinity purification and 
spectrometric analysis (AP-MS) using a combination of pulldown 
assay and direct immunoprecipitation of AMPK-al and -(31 in INS- 
1 beta-cells. After a significance analysis using the SAINT program, 
our interactome comprised 325 and 243 proteins that interacted with 
AMPK-al and -(31, respectively, 196 (51%) of which were novel 
partners. 

Next, we performed a functional classification of our interactome. 
Several proteins that regulated actin organization were linked to 
pancreatic beta-cell functions including, GSIS, beta-cell develop- 
ment, beta-cell differentiation, and cell-cell communication. We also 
validated those interacting proteins that regulated actin cytoskeletal 
organization. Our interactome data constitute a substantial amount 
of new information on AMPK- specific interactions. 

Results 

Overall scheme for profiling of AMPK-al and -01 interactomes. 

To identify the interactomes of AMPK-al and AMPK-(3l in pan- 
creatic beta-cells, we performed a comprehensive AP-MS analysis, 
combining pulldown experiments using tagged recombinant baits 
and direct immunoprecipitation (direct IP) of endogenous target 
proteins (Figure 1). 

Initially, we planned to overexpress 6-myc epitope-tagged AMPK- 
al and AMPK-J31 by transfection in INS-1 cells. However, we could 
not obtain a sufficient amount of myc-tagged proteins to perform 
AP-MS, due to the low transfection efficiency of INS-1 cells 
(Figure 2A). Alternatively, we overexpressed N-terminally 6-myc- 
tagged rat AMPK-al and AMPK-|3l in HEK293T cells and immo- 
bilized them to Dynabeads (Figure 1A, Bait production). Briefly, to 
express the recombinant proteins, full-length rat AMPK-al and -(31 
cDNA was subcloned into modified pcDNA3.1 and transfected into 
HEK293T cells by PEI method 23,24 . 

For the pulldown approach (Figure 1), we used magnetic beads 
(Dynabeads) as the IP material. Recombinant AMPK subunits from 
1 mg of transfected HEK293T cells were immobilized using anti- 
myc-conjugated Dynabeads. To prevent any bias from immobiliza- 
tion of the antibody, anti-myc-conjugated Dynabeads that were 
incubated with mock-transfected HEK293T lysates were used as bait 
in the negative control experiment (Figure 1A). To remove nonspe- 
cific lysate products, INS-1 lysates were precleared with the same 
material as in the IP (Figure IB) 25 . Subsequently, the immobilized 
6-myc-tagged AMPK-al and -(31 were used as bait in pulldown 
assays with 1 mg of INS-1 lysate (Figure 1C). Eluted samples from 



the bait production (Figure 1A) were analyzed by LC-MS/MS to 
confirm the coelution of false-positive interactors and to examine 
homologous interactions in HEK293T cells (Figure ID, HEK293T 
control experiments). 

The proteins that coprecipitated with the myc-tagged baits were 
separated on a 1-DE 10% bis-tris SDS-PAGE gel and subjected to in- 
gel trypsin digestion. Finally, the peptide mixtures were desalted 
using Ci 8 -Stagetips 26 and analyzed on an LTQ Velos linear ion trap 
LC-MS/MS system (Figure IE). Two biological replicates and 3 tech- 
nical replicates were analyzed. 

Further, to increase the coverage and enhance the reliability of the 
AMPK interactome, direct IP was performed without overexpression 
of the target baits (Figure 1). Briefly, native protein complexes were 
isolated directly using antibodies against AMPK-al and AMPK-pT 
from 10 mg of INS-1 cell lysate. Anti-IgG served as the negative 
control. The proteins that coprecipitated with the endogenous baits 
were eluted with acid buffer and subjected to filter-aided sample 
preparation (FASP) 27,28 . Finally, peptide desalting and MS analysis 
were performed as described above (Figure IE). 

Consequently, we generated 4 datasets of the AMPK interactome 
(AMPK-al pulldown, AMPK-al direct IP, AMPK-(3l pulldown, 
and AMPK-[3l direct IP); all experiments were performed with bio- 
logical duplicates and technical triplicates. 

Expression, immobilization, and pulldown were verified by west- 
ern blot (Figure 2). Overexpression of myc-tagged AMPK-al and 
-(31 in HEK293T cells was confirmed by western blot with anti-myc, 
whereas we could not detect any myc-tagged proteins in INS-1 cells 
(Figure 2A). Western blot with AMPK-al and -(31 antibodies was 
also performed to verify the immobilization and coelution of AMPK- 
al and -(31 in the elution fractions of each IP step. AMPK-al and -[31 
were detected in the elution fractions of the target baits but not the 
mock control lane in the pulldown (Figure 2B), indicating that our 
approach reliably identified interactors of the targets of interest. 

Finally, aliquots of the pulldown and immunoprecipitation frac- 
tions were loaded onto an SDS-PAGE gel and stained with silver 
nitrate. Bands that appeared to correspond to AMPK-al and -(31 
and their interacting partners were observed in the target bait IP 
(Supplementary Figures 2A & 2B). 

Identification of novel AMPK-al- and -01- interacting proteins. 

For the pulldown approach, recombinant rat AMPK subunits were 
overexpressed as baits in HEK293T cells, suggesting that some of the 
identified interactors originated from human cells. Thus, the 
proteins that interact with AMPK subunits in INS-1 cells must be 
classified with regard to species. To separate the interacting proteins 
that were derived solely from HEK293T cells, MS/MS spectra were 
searched on the Sorcerer- SEQUEST platform 29 using a target-decoy 
strategy against a concatenated database that contained the human 
IPI database (version 3.87, 91,464 entries), rat IPI database (version 
3.87, 39,925 entries), and their reverse-complements. Conversely, for 
the direct IP approach, MS/MS spectra were searched only in the 
Sorcerer-SEQUEST platform 29 using a target-decoy strategy against 
a concatenated database of the rat IPI database. All search results 
were validated and filtered using Scaffold 4 to establish AP-MS 
datasets at an FDR < 1.0% (Figure 3A). 

Raw data on 36 LC-MS/MS spectra were collected, resulting in 
20,421 unique spectra; 18,377 unique peptides, corresponding to 
1233 unique proteins, were identified at an FDR of 1% 
(Supplementary Table SI). In experiments in which AMPK-al was 
used as bait, an average of 286 and 443 unique proteins were iden- 
tified in the 2 biological replicates of the pulldown and direct IP 
approaches, respectively. With AMPK- (31 as bait, an average of 
179 and 447 unique proteins were identified in the 2 biological repli- 
cates of the pulldown and direct IP approaches, respectively 
(Table 1). 
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Figure 1 | Experimental scheme to examine interactomes of AMPK subunits. The AP-MS procedure comprised 5 steps in series: (A), bait production, 
(B). prey production, (C). pulldown assay, (D). control experiments, and (E). identification. To compare protein interactions, we performed pulldown 
and direct immunoprecipitation (IP). In the pulldown approach, recombinant myc-AMPK was purified from HEK293T lysates using anti-myc- 
conjugated Dynabead baits and linked covalently to matrices for incubation with INS-1 cell lysate (B). The proteins that bind to Myc-tagged AMPK 
proteins were immunoprecipitated and separated on a 1-DE 10% Bis-Tris gel and subjected to in-gel trypsin digestion (C). The digested samples were 
desalted using Ci 8 stage tips. After sample clean-up, the peptides were analyzed by LC-MS/MS. A database search was performed using the Sorcerer 
platform against a concatenated target-decoy database that contained the rat IPI and human IPI databases and their reverse-complements (E). In the 
direct IP, endogenous AMPK subunits were immunoprecipitated in INS-1 cells. Interacting proteins between the control and target AMPK groups were 
analyzed using Scaffold 4. 



Bait recovery of the pulldown approach yielded an average of 294 
and 299 spectral counts with 33% and 57% coverage for tagged 
AMPK-al and -(31, respectively. In addition, bait recovery of the direct 
IP approach generated an average of 42 and 20 spectral counts with 
17% and 11% coverage for tagged AMPK-cxl and -pi, respectively 
(Supplementary Table SI). Because AMPK-al and ~pi form a com- 
plex in cells, we examined whether AMPK subunits were detected in 
our datasets to determine the reliability of our proteomic approaches. 
As shown in Supplementary Figure 3 and Supplementary Table Si, 
nearly all AMPK subunits that were expressed in pancreatic beta-cells, 
such as AMPK-al, AMPK-a2, AMPK-Pl, and AMPK-yl, were iden- 
tified with the target baits, whereas the control baits did not contain 
the spectral counts of AMPK subunits. These results indicate that our 
AP-MS approaches are reliable. 

To remove nonspecific interactors, we performed SAINT ana- 
lysis 30 , as shown in Figure 3A. First, in each experiment, unweighted 
spectral counts of proteins that were identified in the samples that 
corresponded to the control and baits were exported from Scaffold 4. 
Then, to include only high-confidence proteins, proteins had to have 
been identified in all 3 technical replicates with a protein probability 
score s 0.95 to be input into the SAINT analysis (Supplementary 



Table S2). In 4 individual experiments per target bait, identifications 
were accepted for proteins that appeared in at least 2 experiments 
with an average individual probability for SAINT analysis (AvgP 
SAINT) s 0.50. In addition, single-protein hits were considered 
positive at AvgP SAINT > 0.90 (Supplementary Table S3). 

Finally, for the pulldown approach, the amino acid sequences of 
the filtered proteins were examined with regard to species using the 
protein grouping algorithm in Scaffold 4 to sort INS-l-specific inter- 
actors. In Scaffold 4, proteins that are identified with the same 
sequence of a peptide (ie, the sharing peptide) are represented as 
protein groups. Because we used a concatenated human and rat 
database in the pulldown approach, we could verify whether the 
identified proteins were human or rat. For example, if a protein 
group comprises only rat proteins, it clearly matches rat proteins- 
not its human orthologous counterpart. However, if a protein group 
contains a mixture of human and rat proteins, we can conclude that it 
was derived from HEK293T or INS-1 cells. 

Thus, the proteins that were identified in the pulldown approach 
were classified into 3 groups: Rat_only, Rat_Human_Concur- 
rently, and Human_only. To avoid missed interactions and recover 
homologous interactions in the pulldown approach, only proteins in 



SCIENTIFIC REPORTS | 4:4376 | DOI: 10.1038/srep04376 



3 



HEK293T 



INS-1 



o 

o 



Transfection 

73KDa 
41 KDa 



73KDa 
64KDa 

41 KDa 
30KDa 

42KDa 



0. 



CO 

+ 



Q- 



5 

CO 



o 
o 



s 

Q. 



o 



CO 

+ 



C2. 
CL 
< 



CO 

+ 



B 



WB:Anti-Myc 



WB:AMPK-a1 



WB:AMPK-P1 



WB:p-Actin 



AMPK-a1 AMPK-pl 



Transfection - (mock) + 
73KDa 

41 KDa 



• (mock) + 



WB:AMPK-a1 



WB:AMPK-P1 



IP : Anti-Myc 

Figure 2 | Purification of AMPK subunits in control versus target myc-AMPK IP groups. (A) Lysates from HEK293T cells and INS-1 cells expressing 
pcDNA3.1 or Myc-AMPK al/pi were immunoblotted with anti-c-Myc. (B) Lysate from HEK293T cells expressing pcDNA3.1 or Myc-AMPK al/pi 
was incubated with anti-c-Myc-conjugated Dynabeads and covalently coupled to matrices for incubation with INS-1 cell lysate. Bound proteins were 
eluted and subjected to western blot with anti-AMPK-al and -(51. 



the Rat_Human_Concurrently and Human_only groups that had at 
least 99% sequence identity between the rat and human homologs 
were included in the final list by BLAST analysis (Figure 3B). In total, 
325 AMPK-al and 243 AMPK-pT interactors were identified in this 
study and classified as the "INS-1 specific group". In addition, 97 
AMPK-al- and 92 AMPK- (31 -binding partners in the HEK293T 
control experiments (the human species group, Fig. ID) are shown 
in Supplementary Table S4. For further analysis, the HEK293T con- 
trol experiments were processed under the same conditions as with 
the INS-1 -specific group. 

Functional classification of AMPK-specific interactors. To deter- 
mine the functions of proteins that interact with AMPK subunits in 
INS-1 cells, the AMPK-al and -|3l binding partners were grouped by 
Gene Ontology (GO) term using the DAVID bioinformatics 
resource tool into biological process, cellular compartment, and 
molecular function 31 . To identify the GO categories that were 
overrepresented in proteins in INS-1 cells, the INS-1 -specific 
groups and human species group were subjected individually to 
GO analysis (Figure 4). 

By GO analysis, AMPK-al- and -pT -interacting proteins were 
enriched globally in AMPK-mediated processes with similar patterns, 
regardless of species affiliation (Figures 4A and 4B). For example, the 
commonly enriched categories in the 2 groups of the datasets were 
linked to major pathways that are associated with AMPK, such as 
metabolism 3,5 , translation 32 , and biosynthesis. Also, most proteins 
were distributed broadly throughout the cytosol, mitochondria, nuc- 
leolus, and ribosome. AMPK-al- and -[51 -interacting proteins that 



mediate the organization and regulation of the actin cytoskeleton 
were significantly enriched in the INS-l-specific versus human species 
group (Figures 4A and 4B). In addition, many proteins in the INS-l- 
specific group were involved in transport, the cell cycle, and cyto- 
skeletal protein binding. The similarity in patterns between AMPK-al 
and -[31 interactors in the GO analysis might be attributed to the 
structural conformation of the heterotrimeric AMPK complex, which 
contains AMPK-al and -(31. All information on the GO annotations 
is listed in Supplementary Table S5. 

To examine the signaling pathways in our AMPK interactome, we 
analyzed our data using the KEGG pathway database (http://www. 
genome.jp/kegg/) and DAVID bioinformatics tool 31 . A total of 381 
AMPK-al- and -(31 -interacting proteins in the INS-l-specific group 
fell into the following KEGG pathway categories: Ribosome, 
Spliceosome, Huntington disease, Parkinson disease, Adherens junc- 
tion, and Regulation of actin cytoskeleton (Supplementary Table S6). 
Notably, KEGG pathway analysis of the human species groups 
demonstrated that only proteins in the INS-l-specific groups were 
significantly enriched with regard to regulation of the actin cytoske- 
leton (Supplementary Table S6). 

Fifteen proteins regulated organization of the actin cytoskeleton, 
which is consistent with the GO analysis. In addition, based on a 
literature search 33 , Myll2b was linked to actin cytoskeletal organiza- 
tion. Consequently, 16 proteins that were associated with actin cyto- 
skeletal organization were obtained, based on the functional 
classification (Table 2). 

In pancreatic beta-cells, the actin cytoskeleton has significant 
functions in many processes, including GSIS, cell-cell adhesion, 
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and proliferation 22,34,35 . Although whether AMPK interacts directly 
with components in the actin cytoskeleton is unknown, recent stud- 
ies have indicated that AMPK regulates its organization to mediate 
cell mobility and migration 36 3S . Based on the novelty of these results, 
we focused on the interactions of proteins that are associated with 
actin cytoskeletal organization. 

Interaction of proteins related to actin cytoskeletal organization 
with AMPK. To examine the interactions between proteins that 
regulate actin cytoskeletal organization and AMPK, we generated a 
protein interaction network of 140 proteins from the KEGG pathway 
database using STRING 9.1 39 and Cytoscape 40 . The evidence-based 
global STRING-generated protein network 39 showed high connec- 
tivity between actin cytoskeleton-regulating proteins (Figure 5A). 
However, the relationship between AMPK and actin cytoskeletal 
organization was not well defined in the protein interaction 
network. Thus, we added 28 direct interactions between 17 actin 
cytoskeleton-related proteins and AMPK-al/-(il to the protein 
interaction network (Figure 5A & Table 2). 

Of the interactions that were detected by AP-MS, the binding of 
Iqgapl, Gsn, and Vim to AMPK-al and -(31 was confirmed by 
coimmunoprecipitation using 6-myc-tagged AMPK-al and -(31. 
As shown in Figure 5B, IQGAP1 and gelsolin (Gsn) coimmunopre- 
cipitated with both AMPK subunits. Vimentin (Vim) was used as the 
positive control, based on the robust interaction between vimentin 
and AMPK. 

The interaction with Myh9, IQGAP1, Rhoa, and Racl was vali- 
dated by direct coimmunoprecipitation in INS-1 cells (Figure 5C). 
The interaction between AMPK-al and -(31 was tested as the positive 
control. In total, we validated the interactions of 6 proteins (Gsn, 
Iqgapl, Vim, Myh9, Rhoa, and Racl) that were linked to regulation 
of the actin cytoskeleton, demonstrating that AMPK-al and -(31 
bind to the actin cytoskeleton and interact physically with regulators 
of the actin cytoskeleton organization. 



Discussion 

The aim of this study was to identify AMPK-interacting proteins in 
pancreatic beta-cells. For the AP-MS analysis, we first examined the 
expression of tagged baits in HEK293T cells, a cell type that has high 
transfection rates, because INS-1 beta-cells have low transfection effi- 
ciency of DNA vectors. Although we tried to exclude binding partners 
that arose solely from the HEK293T cells through a negative control 
experiment and species affiliation analysis with a combination data- 
base (human and rat) search, it is possible that false-positive inter- 
actors and overlooked interactions existed in our pulldown approach. 

To verify the reliability of the species affiliation analysis in the 
pulldown approach, we performed several preliminary control 
experiments. In generating the tagged AMPK baits (Figure 1A), 
the proteins that were considered to have originated from 
HEK293T cells were analyzed by LC-MS/MS without incubation 
with INS-1 cell extracts (Figure ID), whereas the baits was incubated 
with INS-1 cell extracts in the pulldown assay (Figures IB & 1C). 
These experiments indicated that that species affiliation analysis 
distinguished the origins of proteins in the pull-down of INS-1 cell 
extracts with myc-AMPK-al and -(31 (Supplementary Figure 4A). 

Further, BLAST analysis was performed to prevent missed inter- 
actions and increase the sensitivity in the pulldown assay. Only pro- 
teins with at least 99% sequence identity between the human and rat 
orthologs in the Rat_Human_Concurrently and Human_only 
groups were included in the final list. To ensure the specificity of 
the BLAST analysis, vimentin (98% sequence identity between 
human and rat), an interactor of AMPK-a and -(3 21,41 , was used as 
the positive control. Despite its sequence identity (98%), we con- 
firmed a homologous interaction between vimentin and AMPK sub- 
units (Figure 5B), suggesting that the BLAST analysis could rescue a 
homologous interaction in the pulldown approach. 

As shown in Supplementary Figure 4B, our BLAST analysis 
recovered approximately 50% of proteins that were identified as 



SCIENTIFIC REPORTS | 4:4376 | DOI: 10.1038/srep04376 



5 



o- 
o 
d 

Al 



< 



D 

"5 



0) 



0) 

to 



D 



0 
Q_ 



O 
D 



_0 
D 



3? E 2 
Z 



J5 



i § 

?,<] 



— >- 



h°e a. 



CD </> 

P-| 



d^ 3 



2 i> 
z 



P| 



5 D 



^z£ 



p-5 



6^ 



— cn co CN co co cn 

KOKOKO^O 

"coio^t^tO'OO 
cn ' 1 ' 1 1 1 ' 



CO CN t\ CO 
IO O CO CO 



-St "St — — 

O K CO CO 



COCOU")-^-CNOK~0 
KOKOiOO"^co 
(N ■— ■— 1— 



KCOOKCOK-^O 
Kco-^-LOKKO — 
cocooo'Onoo 



CNCNKCOOOOK 
OCNCMO^tCN^t^ 
OCOiOCOCNOCO-^- 
CN ■ — < — sf lOM CN 



cocsst-m^fcno 

CNCNOOOOKK 



CNKOOOOCNCO 
COOOKCNOOK 

CN ■ — ' — CO ^t CN CN 



-^■COCNOK-^-^O — 
O^tCN^tK^tCO^t 

— — CO K CO » — O 

(N •— — CO CO CN — 



OCOOOI^t-K 
OO^HOCOCO-OiO 



CNCOOT>OK->KCO 
lOCNCMCOO^tCO^t 
COCNCM^iOOCNCN 



COCNCN^iOOCNCN 



oo--^o.n-ocN 
KO-O-COO-JtM 
CNCN^^CNCN^^ 



&&&&&&&& 
<<<<<<<< 



i- CN 

O O 

IS IS 

1 1 

a. a. 



<N 

0 
bq 
c 

0 

TJ 

0 — — 



0 

3 
c 

0 

TJ 



3 

a. 



3 



Human_only and Rat_human_concurrently after the filtration per 
the 99% sequence identity criterion. Moreover, the similar distri- 
bution of sequence identities in the HEK293T control experiments 
demonstrates that the BLAST analysis was not biased against the 
pulldown assay (Supplementary Figure 4B). 

To overcome the limitations of the pulldown assay and increase 
coverage of the AMPK interactome with high confidence, direct IP 
was performed using antibodies against AMPK-al and -|3l in INS-1 
cells without overexpression of baits. In addition, the pulldown and 
direct IP experiments were repeated independently with each bait 
and processed using the SAINT algorithm 30 . The SAINT statistical 
platform allowed us to assign confidence scores to interacting pro- 
teins from our AP-MS experiments. SAINT uses spectral counts of 
identified proteins to derive a probability score of a true interaction 
and discerns true interaction partners from background noise 30 . 

Thus, we used a stringent AP-MS strategy that was complemented 
by bioinformatics, SAINT analysis, and affinity purification-com- 
bining pulldown assay and direct immunoprecipitation-to draft the 
AMPK-al and -(31 interactome in INS-1 pancreatic beta-cells. 
Because we combined the Human and Rat IPI databases in the pull- 
down assay, the IPI accession numbers of the identified proteins were 
converted to gene names (symbols) and subjected to comparison 
analysis. 

The Venn diagrams in Supplementary Figures 5A show the over- 
lap in unique proteins between the pulldown assay and direct IP for 
AMPK-al and -f3l. These results demonstrate that our approaches 
provide complementary coverage and that their combination effects 
comprehensive coverage of the AMPK-al and -(31 interactomes. 
Although the overlap between datasets was low due to biological 
and technical variations between the approaches and the use of 
low- resolution mass spectrometry 42 , 201 (62%) and 116 (48%) pro- 
teins were identified in at least 2 of 4 biological replicates for AMPK- 
al and -pi, respectively. Notably, between the AMPK-al and 
AMPK- [31 datasets that were merged from the approaches, there 
was 58% to 77% overlap (187 of 325 and 243 proteins in AMPK- 
al and AMPK-|3l each, respectively), which is better compared with 
the individual methods (Supplementary Figure 5B). Although the 
overlap between experiments was not high, the approaches enhanced 
coverage and identified distinct segments of the AMPK interactome. 

We compared our AMPK-al and -(31 interactome with existing 
data. Although several partners of AMPK-al and -[31 have been 
reported 21,41,43 , no global interactome analysis has focused on 
AMPK subunits. In particular, none of the well-characterized part- 
ners of these subunits have been identified in pancreatic beta-cells, 
based on mass spectrometry. 

Large-scale proteomic data on AMPK interactors have been 
reported by Behrends et al. 21 , Pilot-Storch et al. 43 , and Varjosalo 
et al. 41 in human cell lines. The 614 proteins in these reports were 
compared directly with our AMPK-al and -(31 interactors, based on 
gene name (symbol) (Supplementary Figure 6A). Of the 381 proteins 
in our list, there were 145 proteins with the same gene name in the 3 
datasets. In addition, we compared our interactome with a public 
database 44 . A total of 253 interactions for AMPK-al and - (31 were 
obtained from the Protein Interaction Network Analysis (PINA) 2.0 
database 44 , 41 of which were common between our interactome and 
the database (Supplementary Figure 6B). Further, our data were 
compared with the results of the protein-protein interaction predic- 
tion by the PrePPI server 45 , which included 736 and 675 interactions 
of AMPK - al and -(31 in human cells, respectively, 63 of which were 
found in our data (Supplementary Figure 6C). Collectively, 185 
(48%) proteins in our interaction data were identified in other stud- 
ies, and 196 (52%) proteins were postulated to be novel 
(Supplementary Figure 6D). All previously reported proteins that 
associate with AMPK are summarized in Supplementary Table S7. 
Considering all comparison analyses, our dataset contains a substan- 
tial amount of new information on AMPK-specific interaction. 
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Figure 4 | Functional GO annotation of AMPK-al and AMPK-pl interactome. All identified proteins-325 AMPK-al- and 243 AMPK- pi -interacting 
proteins-were grouped, based on GO terms. GO terms enriched from proteins identified in INS-1 specific group were classified as specific, whereas 
those from proteins in the INS-1 -specific and Human species groups were classified as common. Fold-enrichment for a particular GO term describes the 
ratio between the numbers of genes belonging to a specific GO term and the total number of genes in the genome with at least 1 GO annotation. 
Enriched GO terms for biological process (GOBP), cellular component (GOCC), and molecular function (GOMF) are shown. In AMPK-al (A) and 
AMPK-|3l (B), GO terms associated with the regulation of actin cytoskeletal organization were significantly enriched. 



Because our baits are the catalytic and regulatory subunits of the 
AMPK complex, we expected some of its partners to be substrates for 
it. Thus, we determined whether the identified proteins were phos- 
phorylated in tissues or cells and if their phosphorylation was 
AMPK-dependent. Our primary sources for this step were large 
phosphoproteome repositories, based on mass spectrometry, includ- 
ing the PhosphoSitePlus database 46 , PHOSIDA database 47 , and 



several phosphoproteomic studies in human, mouse, and rat 48-53 . 
Moreover, we predicted AMPK-dependent phosphorylation sites 
using GPS2.1 54 to examine the connections between the identified 
binding partners and AMPK. 

As shown in Supplementary Figure 7A and Supplementary Table 
S7, more than 92% of identified proteins were phosphorylated, over 
78% of which (299 proteins) were predicted to be AMPK-dependently 



Table 2 Novel interactors of AMPK associated with regulation of acti 
Gene Symbol Protein Name 
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Molecular Weight 


Species 


baits 


Novel 


Iqgapl 


IQ motif containing GTPase activating protein 1 


1 89 kDa 


Rattus norvegicus 


al/pl 
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Arpc4 


actin related protein 2/3 complex, subunit4 


20 kDa 


Rattus norvegicus 


al 
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Actgl 


actin, gamma 1 


42 kDa 


Rattus norvegicus 


al/pl 
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Actn4 


actinin alpha 4 


1 05 kDa 


Rattus norvegicus 


al/pl 
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Actnl 


actinin alpha 1 


1 03 kDa 


Rattus norvegicus 


al/pl 
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Cfll 


cofilin 1 , non-muscle 


19 kDa 


Rattus norvegicus 


al 
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Cfl2 


cofilin 2, muscle 


19 kDa 


Rattus norvegicus 


PI 
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Gsn 


gelsolin 


86 kDa 


Rattus norvegicus 


al/pl 
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MyhlO 


myosin, heavy chain 10, non-muscle 


234 kDa 


Rattus norvegicus 
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Myhl4 


myosin, heavy chain 14 


232 kDa 


Rattus norvegicus 


PI 


X 


Myh9 


myosin, heavy chain 9, non-muscle 


226 kDa 


Rattus norvegicus 


al/pl 
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Pppl ca 


protein phosphatase 1 , catalytic subunit, alpha isoform 


38 kDa 


Rattus norvegicus 


al/pl 
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Ppplcb 


protein phosphatase 1 , catalytic subunit, beta isoform 


37 kDa 


Rattus norvegicus 


PI 
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ras homolog gene family, member A 


22 kDa 


Rattus norvegicus 


al 


O 


Racl 


ras-related C3 botulinum toxin substrate 1 


21 kDa 


Rattus norvegicus 


al/pi 


O 


Myll2b 


myosin regulatory light chain MRLC2 


20 kDa 


Rattus norvegicus 


al/pl 
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Vim 


Vimentin 


54 kDa 


Homo sapiens 


control 


X 
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Figure 5 | Interaction network for proteins associated with AMPK complex in regulation of actin cytoskeleton. (A) Network analysis of AMPK 
interactome, showing that the AMPK complex interacts directly with multiple proteins associated with the regulation of the actin cytoskeleton. Proteins 
and protein-protein interactions associated with regulation of the actin cytoskeleton were extracted from the KEGG pathway and STRING 9. 1 databases. 
All input proteins associated with actin regulation are depicted as small gray spheres. Interacting proteins in our dataset are depicted as green spheres, and 
interacting proteins validated by Co-IP are shown as yellow spheres. Bait proteins are shown as red hexagons. Protein-protein interactions extracted from 
the STRING database are shown as grey lines. Interactions identified in this study are shown as green lines. Co-IP-validated interactions are shown as red 
dotted lines. (B) Validation of protein-protein interactions by pulldown approach. A novel interacting protein (Iqgapl) was validated by co-IP. Well- 
known proteins (Gsn and Vim) were also validated. Vim was used as a positive control to confirm homologous interactions. (C) Validation of protein- 
protein interactions using direct IP. AMPK-al, AMPK-fil, Myh9, Iqgapl, Rhoa, andRacl were validated by co-IP. The co-IP results for novel interactors 
(Myh9, Iqgapl, Rhoa, and Racl) that were associated with regulation of the actin cytoskeleton corresponded well with the AP-MS results. 



phosphorylated, using GPS2.1 with moderate confidence (Supple- 
mentary Figure 7B). Further, of 1410 phosphorylation sites in the 
299 proteins, 454 sites in 169 proteins were identified in earlier MS- 
based phosphoproteome studies (Supplementary Figure 7C and 
Supplementary Table S8). Detailed information on the prediction 
analysis using GPS2.1 is shown in Supplementary Figure 7D. For 
AMPK -al and -fil, nearly 80% of interactors were predicted to have 
phosphorylation sites with medium confidence (Supplementary 
Figure 7D). 

Considering proteins that have been reported to associate with or 
be linked functionally to AMPK, our comparison analysis suggests 
that 92% (349/381) of identified proteins are related to AMPK, cor- 
roborating the functional correlation of these proteins with AMPK 
and validating the reliability of our approach for discovering novel 
AMPK interaction partners. 

We screened the AMPK- a 1 and -(31 interactomes through 
functional grouping, based on their affiliation, and sorted AMPK 
interactors that were derived from pancreatic beta-cells. Most 
AMPK-al- and -(3l-binding proteins in the INS- 1 -specific group 
regulated actin cytoskeletal organization. Because a principal func- 
tion of pancreatic beta-cells is to secrete hormones on stimulation, 
such as insulin through actin cytoskeleton remodeling 55-56 , the pro- 
teins that are associated with the stimulus could be upregulated or 
activated in beta-cells versus other cell types, which might explain the 
GO analysis results. 



Our bioinformatics analysis and literature search revealed that 16 
proteins are linked to the regulation of actin cytoskeletal organiza- 
tion (Table 2)-Myh9, MyhlO, Myhl4, and Myll2b are associated 
with myosin class II (myosin II) activation and, as actin-binding 
proteins, are thus central in cell adhesion, migration, and division 57 . 
Myh9, MyhlO, and Myhl4 are 3 members of nonmuscle myosin II 
(NMII), including NMIIA, NMIIB, and NMIIC. Whereas MyhlO 
and Myhl4 form complexes with AMPK in human cells 41,58 , our 
AP-MS and co-IP data constitute novel findings of a link between 
AMPK and Myh9. Recent studies have suggested that NMII mod- 
ulates GSIS through actin and focal adhesion remodeling in pancre- 
atic beta-cells 59 . Further, AMPK regulates cell morphology and 
secretion by phosphorylating Myhll2b (known as MRLC) 19,60 . 
Although AMPK controls myosin II directly 61,62 and indirectly 19,63 
in pancreatic beta-cells, role of AMPK in regulation of myosin 
remain unresolved. Our data support a model in which AMPK med- 
iates myosin II-associated beta-cell functions, including the regu- 
lation of GSIS. 

Notably, several proteins (IQGAP1, gelsolin, cofilin, Racl, and 
RhoA) are related to reorganization and remodeling of filamentous 
actin (F-actin) 56 , which regulates insulin secretion in beta-cells 56 . As 
scaffolding proteins and F-actin-binding protiens, IQGAP1 and gel- 
solin regulate cytoskeletal reorganization 64 and are also linked to 
pancreatic beta-cell functions, such as insulin secretion 55,65,66 , beta- 
cell development 67,68 , and cell-cell communication 66,69 . Whereas gel- 
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solin interacts with AMPK-p in human cells 21 , we report a novel 
interaction between AMPK and IQGAP1. Both proteins were con- 
firmed to interact with AMPK-al and -pi by co-IP and western blot 
(Figure 5B). 

Cofilin-1 and Cofilin-2 are associated with glucose homeostasis, 
regulating insulin-stimulated GLUT4 translocation in muscle cells 70 
and pancreatic beta-cells 71 . Racl is a member of the Rho family of 
GTPases, which govern cytoskeletal organization 72 . Recent evidence 
also implicates an association between AMPK and Racl in insulin 
secretion by pancreatic beta-cells 56 ' 72 . RhoA is a small GTPase that 
regulates the actin cytoskeleton during the formation of stress fibers 
and is phosphorylated by AMPK 73 . Interestingly, RhoA controls sev- 
eral beta cell functions including beta-cell spreading and actin cytos- 
keleton dynamics through the Rho-ROCK pathway 74 . 

Because most of the direct interactions between AMPK and reg- 
ulators of actin reorganization have not been reported in pancreatic 
beta-cells, our interactome is a resource that can be used to examine 
novel functions AMPK in regulating actin remodeling in pancreatic 
beta-cells. Further, the validation of 5 proteins (Gsn, Iqgapl, Myh9, 
Rhoa, and Racl) should help us determine the function of AMPK- 
actin cytoskeleton assembly in regulating GSIS. 

In conclusion, we have identified many putative partners of 
AMPK-otl and -pi by combining affinity chromatography and mass 
spectrometry. With our stringent filtering criteria, our interactome 
comprised 325 and 243 proteins that interacted with AMPK-al and - 
pi, respectively. Our study confirms 185 previously described pro- 
teins and identified 196 novel proteins that bind AMPK in our 
AMPK interactome datasets. Our data reinforce a model in which 
AMPK governs many functions that are linked to the regulation of 
actin organization. Specifically, several such proteins are closely 
associated with pancreatic beta-cell functions, including GSIS, 
beta-cell development, beta-cell differentiation, and cell-cell com- 
munication. Consequently, our interactome can be used to perform 
more detailed biochemical analyses of putative substrates and effec- 
tors of AMPK. 

Methods 

Reagents and materials. HPLC-grade acetonitrile (ACN), HPLC-grade water, 
HPLC-grade methanol (MeOH), hydrochloric acid (HQ), and sodium chloride 
(NaCl) were obtained from DUKSAN (Kyungkido, Korea). Acetic acid was 
purchased from TEDIA (Fairfield, OH, USA). Brilliant Blue G-250 and a Bio-Rad 
protein assay kit were purchased from Bio-Rad (Hercules, CA, USA), and Complete 
Protease inhibitor cocktail was purchased from Roche (Mannheim, Germany). Other 
reagents, including 2-mercaptoethanol (|3-ME), ammonium bicarbonate (ABC), 
ammonium persulfate (APS), EDTA, iodoacetamide (IAA), magnesium chloride, 
MES hydrate, potassium chloride, sodium bisulfite, sodium deoxycholate (SDC), 
sodium pyrophosphate, and trifluoroacetic acid (TFA) were purchased from Sigma- 
Aldrich (St. Louis, MO, USA). 

Cell culture. The INS-1 rat insulinoma cell line was grown in RPMI 1640 medium, 
containing 10% heat- inactivated fetal bovine serum, 100 units/ml antibiotics, 1 mM 
sodium pyruvate, 10 mM HEPES (pH 7.4), and 50 uM (J-ME, at 37 C in a 5% C0 2 
incubator. 

HEK293T cells were maintained in Dulbecco's modified Eagle's medium (DMEM) 
with 10% heat- inactivated fetal bovine serum and antibiotics (100 units/ml penicillin 
and 100 mg/ml streptomycin). 

Plasmid construction and transient transfection. The full-length coding regions of 
rattus AMPK-al and AMPK-pl were generated by PCR of first-stand cDNA from 
total RNA of INS-1 cells. Each construct was cloned into the eukaryotic pcDNA3.1-6- 
myc-tagged expression vector (Supplementary Figure 1A). Information for 
constructs and primers is shown in Supplementary Figure IB. Cloned plasmids were 
confirmed by sequencing. Plasmids were prepared using the Plasmid Midi kit per the 
manufacturer's protocol. 

The plasmids were transfected into HEK293T cells by PEI method 23,24 . Prior to 
transfection, 12 ug of plasmid in 2 ml serum-free media were mixed with 36 ul PEI 
solution (1 ug/ul) and incubated for 15 min at room temperature. The solution was 
added to 1 X 10 7 HEK293T cells in a 150-cm 2 dish. Cells were incubated for 4 hours at 
37°C in a C0 2 incubator. After being washed with 10 ml of serum-free media twice to 
remove excess PEI, the cells were incubated with 20 ml of serum-containing media 
overnight at 37°C. Transiently transfected cells were washed 3 times with ice-cold 
PBS and pelleted. The pellets were snap-frozen in liquid nitrogen and stored at — 80"C 
until further analysis. 



Preparation of cell lysates. HEK293T and INS-1 cells were lysed in 200 ul modified 
RIPA buffer, comprising 150 mM NaCl, 50 mM Tris-Cl pH 7.4, l mM EDTA, IX 
protease inhibitor cocktail, 0.1 mM PMSF, 0.1% SDC, and 1% NP-40. Cells were 
disrupted by sonication and centrifuged at 15,000 rpm for 40 min at 4°C to remove 
cell debris. Protein concentration was measured by BCA assay. Anti-myc western blot 
analysis was used to monitor recombinant myc-AMPK expression. 

Pulldown assay using anti-myc coupled to magnetic beads. Prior to 
immunopurification, 1 mg each of HEK293T and INS-1 cell lysate was precleared 
using protein G Dynabeads at 200 : 1 (cell lysate to beads) for 4 hr at 4°C on a rotary 
device. First, 50 ul of protein G Dynabeads was washed with PBS and incubated for 
1 h with 5 ug anti-myc at room temperature. After washes with PBS, the recombinant 
6-myc-tagged AMPK subunits were immobilized from precleared HEK293T cell 
lysates onto anti-myc -conjugated Dynabeads at 4°C for 4 h. 

For the pulldown, the beads were washed 3 times with 500 ul RIPA buffer and 
mixed with precleared INS-1 cell lysate overnight at 4°C, After nonspecific proteins 
were washed away, the immunopurified mixtures were incubated with elution buffer 
(0.1 M citrate pH 2.5) at room temperature for 5 min. In the control experiment, the 
same immunopurification was performed using lysates from mock-transfected 
HEK293T cells in the immobilization step. 

SDS-PAGE separation, in-gel digestion, and desalting. Immunopurified fractions 
of each IP were resolved in SDS sample buffer and denatured at 95 C for 5 min. After 
protein separation by 10% Bis-Tris SDS-PAGE and staining with Brilliant Blue G- 
250, each lane of the gel was excised, cut into 6 slices, and incubated with 500 ul 
200 mM ABC and 50% ACN for destaining. The gel slices were then dehydrated and 
rehydrated with 200 ul 100% ACN and 200 ul of 0.1 M ABC, respectively. Disulfide 
bonds were reduced with 200 ul 10 mM DTT. Then, the proteins were alkylated with 
50 mM IAA for 30 min at RT in the dark. 

The gel pieces were dried in a speed vacuum and incubated overnight at 37°C with 
sequencing- grade modified trypsin (Promega, Fitchburg, WI, USA) at an enzyme-to- 
protein ratio of 1 : 100 (w/w). The resulting peptides were extracted sequentially from 
the gel slices with 100 ul40%ACN/50 mMABCandlOO ul 80% ACN/0.1% TFA by 
sonication for 15 min at each stage. After all supernatants were combined, the pep- 
tides were dried in a speed vacuum. 

The dried peptide mixtures were dissolved in 100 ul 0.1% TFA and desalted using 
homemade C 18 -StageTips as described 26,27 . Briefly, C 18 -StageTips were prepared by 
packing POROS 20 R2 material (Applied Biosystems, Foster City, CA, USA) into 200- 
ul yellow tips on top of C 18 Empore disk membranes. The C 18 -StageTips were washed 
with 100 ul 100% ACN and equilibrated 3 times with 100 ul 0.1% TFA using a 
syringe. After the samples were loaded, the C 18 -StageTips were washed 3 times with 
0.1% TFA and eluted with 100 ul of a gradient of elution buffers, containing 0.1% 
TFA and 40%, 60%, or 80% ACN. All eluates were combined, dried on a speed 
vacuum, and stored at — 80°C until further analysis. 

Direct immunoprecipitation (IP). To identify interactors of AMPK-al and AMPK- 
pi, we performed direct immunoprecipitation assay using INS-1 cell lysate. After 
protein G Dynabeads were coated with rabbit monoclonal AMPK-al or AMPK-|3l 
antibody as bait or normal rabbit IgG as negative control, the beads were mixed with 
10 mg of precleared INS- 1 cell lysate. The immunopurified mixtures were eluted with 
0.1 M citrate buffer. One-third of the volumes of the eluate was digested by FASP as 
described 27 . The resulting peptides were desalted as described above. The remaining 2 
volumes were separated by 1-D electrophoresis on 8-15% SDS-PAGE gels and 
immunoblotted with antibodies of novel interacting partners for validation. 

LC-MS/MS analysis. Peptide mixtures were analyzed on an EASY nano LC (Proxeon, 
Odense, Denmark), interfaced with an LTQ Velos mass spectrometer (Thermo 
Electron Corporation, San Jose, CA, USA), as described 27 with some modifications, 
nano LC was operated in the 2-column system with a trap column (75 um I.D. X 
4 cm) and analytical column (75 um I.D. X 15 cm) that were packed in-house with 
Cis resin (Magic Ci 8 -AQ, 5 um, 100 A). Solvent A was 0.1% formic acid and 2% 
ACN, and solvent B was 0.1% formic acid and 98% ACN. Fifty milliliters of each 
sample, dissolved in 50 ul solvent A, was injected into the trap column at 5 ul/min. 
Peptides were eluted with a gradient of 2% to 40% solvent B over 85 min, followed by 
a gradient of 40% to 90% for 15 min and 90% over 5 min at 500 nl/min. 

The ion spray voltage was set to 1.8 kV in the positive ion mode, and the tem- 
perature of the heated capillary was 320°C. A cycle of 1 precursor MS survey spectrum 
(300-2000 m/z) was acquired in the profile mode. MS/MS scans were taken in the 
linear trap in a data-dependent manner for the 10 most abundant signal precursor MS 
ions. All CID MS/MS spectra were acquired using the following parameters: 35% 
normalized collision energy; ion selection threshold of 500 counts; activation Q of 
0.25; and activation time of 30 ms. Dynamic exclusion was performed with a repeat 
count of 1, 30-s repeat duration, exclusion list size of 50, exclusion duration of 60 s, 
and ±1.5 m/z exclusion mass width. Overall, 3 technical replicates were analyzed for 
each dataset. 

Data processing. Raw MS files were processed using the Sorcerer-SEQUEST 
platform 29 as described with some modifications 27 . MS/MS spectra was examined 
using a target-decoy database search strategy against a concatenated forward/ 
reversed version of the International Protein Index (IPI) rat (v 3.87, 39,925 entries) 
and human (v 3.87, 91,464 entries) databases, supplemented with the protein 
sequences of bovine serum albumin, trypsin, and 6-myc-tagged AMPK-al and 
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AMPK- [31. The database search parameters were: full enzyme digest using trypsin 
(After KR/-) with up to 1 missed tryptic cleavage; a parent ion mass tolerance of 
2.0 Da {average mass); a fragment ion mass tolerance of 0.8 Da (monoisotopic mass); 
fixed modification of 57.02 Da on cysteine for carboxyamidomethylation; and 
variable modifications of 15.99 Da on methionine for oxidation. 

The search data for each experimental scheme were merged and validated using 
Scaffold 4. Collected datasets were filtered, based on peptide probability, protein 
probability, and SEQUEST scores, to achieve an estimated false discovery rate (FDR) 
< 1.0%. The filter criteria were as follows: for charge states of 1 + , 2+, and 3 + , Xcorr 
scores should be greater than 1.5, 2.5, and 3.5, respectively; deltaCN should be larger 
than 0.1; peptide probability and protein probability should exceed 0.95; and the 
minimum number of unique peptides was set to 2. 

To identify proteins that interacted specifically with AMPK subunits and eliminate 
false interactions from the negative controls, significance analysis of interactome 
(SAINT) was performed, based on spectral counts, as described 30,75,76 . The probability 
scores of the bait and prey proteins were calculated as the average of the probabilities 
in individual replicates (AvgP). Proteins with AvgP > 0.9 in 1 biological replicate or 
that were detected in at least 2 of 4 biological replicates with AvgP > 0.5 were likely 
interactors 75,76 . 

Bioinformatics analysis. Gene ontology was performed using DAVID 
Bioinformatics Resource 6.7 31 . Interactors were categorized into cellular 
compartment, molecular function, and biological process terms. Only enriched GO 
terms with a p- value < 0.05 were selected. Pathway analysis was performed using the 
KEGG (Kyoto Encyclopedia of Genes and Genomes) Pathways database (http:// 
www.genome.jp/kegg/). A protein interaction network was generated using STRING 
9.1 39 and visualized with CytoScape 40 . AMPK-specific phosphorylation sites were 
predicted using GPS2.1 54 . The prediction threshold values for GPS2.1 were set to high 
(3.081) and medium (2.324). Based on the prediction values, phosphorylation sites 
were classified as high confidence and medium confidence. Phosphorylation sites 
with low prediction values constituted the no-prediction group. 

Immunoblot analysis and antibodies. Whole HEK293T and INS-1 cell lysates were 
prepared from 1 X 10 7 cells in 200 ul RIPA buffer. Pulldown assay and direct IP assay 
were performed using Dynabeads as described above. Control experiments were 
performed using lysates from mock-transfectedHEK293T cells. 

The mixtures were incubated overnight at 4° C with rotation to capture the immune 
complex. The beads were pelleted by brief centrifugation and washed 3 times with 
RIPA buffer. Bound proteins were eluted with 100 ul elution buffer, consisting of 

0. 1 M citrate. The resultant proteins were separated on an 8% to 15% SDS-PAGE gel 
and immunoblotted with antibodies against AMPK-al and AMPK-J31 (Cell Signal, 
Boston, MA, USA), MYH9 (sc-47199), IQGAP1 (sc-10792), gelsolin (sc-48769) and 
RhoA (sc-179) (Santa Cruz Biotechnology, Santa Cruz, CA, USA), and Racl (05-389) 
(Millipore, Billerica, MA 01821, USA). 
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